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We study statistical properties of a class of band random matrices which naturally appears in 
systems of interacting particles. The local spectral density is shown to follow the Breit-Wigner 
distribution in both localized and delocalized regimes with width independent on the band/system 
size. We analyse the implications of this distribution to the inverse participation ratio, level spacing 
statistics and the problem of two interacting particles in a random potential. 
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Intensive investigations of band random matrices 
(BRM) have been done during last years Differ- 
ent regimes corresponding to localized and delocalized 
wave-functions have been studied numerically and ana- 
lytically and it has been shown that the transition from 
one regime to another can be described by one scaling 
function depending on the ratio of the localization length 
in the infinite system I ~ b 2 to the size of the matrix N, 
where the parameter b determines the size of the band 
2b + 1. Similar types of matrices appear in such physical 
systems as quasi one-dimensional disordered wires and 
such models of quantum chaos like the kicked rotator 
that gives additional grounds for investigation of BRM. 

The above BRM can be also considered as a reason- 
able model of one-particle localization in a disordered 
wire of finite size p],0]. However recent investigations 
of two interacting particles (TIP) in a random potential 
H showed that another type of BRM naturally appears 
in interacting systems. Indeed for interacting particles 
there is one preferential basis which corresponds to eigen- 
states without interaction. In this basis, the total Hamil- 
tonian is the sum of a diagonal matrix, with elements 
given by the sum of one-particle energies, and a BRM, 
which describes interaction induced transitions between 
eigenstates of the non-interacting problem. The first in- 
vestigations of such superimposed BRM (SBRM) allowed 
to find the dependence of the localization length l s b on 
the amplitude Wb of large fluctuations on the diagonal 
and to obtain the localization length l c for two-particles 
coherent propagation in a random potential on a distance 
much larger than one-particle localization length l\ ||. 

While from the TIP model it is clear that matrices 
with preferential basis should describe interesting physi- 
cal effects in interacting systems, only few investigations 
in this direction have been done up to now [QJ^]. In this 
paper we investigate the properties of such matrices in 
particular the local spectral density and the inverse par- 



ticipation ratio (IPR). Due to the close connection be- 
tween the SBRM and the TIP problem, the obtained 
results can also be used for the latter case. 

The matrix we study is the sum of a random diagonal 
matrix and a conventional BRM : 



Hn,n' = VnS n ,n> + Cn,n' / V^b + 1 



(1) 



with -W b < Vn < W b , -1 < Cn,n> < 1 for \n ~n'\<b 
and £ n ,n' = elsewhere. The connection with the TIP is 
given by b ~ If and W b ~ A\/l\V/U in terms of the inter- 
action strength U and the one-particle energy bandwidth 
4V ', l\ 3> 1. This matrix describes a one-dimensional 
two-particle Anderson model, with on-site interaction U, 
in the basis of non-interacting eigenstates. In |^] it was 
shown that the eigenstates of (1) are localized with lo- 
calization length l sb « b 2 /2W b 2 for 1 < W 6 < Vb. This 
leads to an enhancement of the length of coherent TIP 
propagation l c — l s b/h ~ l 2 (U/V) 2 /32 independent on 
the sign of interaction. 

Our numerical investigations of SBRM (1) show that, 
in addition to the standard exponentially localized form, 
the eigenstates are also characterized by large amplitude 
fluctuations of probability on nearby sites. A typical ex- 
ample of such an eigenstate is presented in Fig.l. The 
spike eigenstate structure is clearly noticeable. 

This implies that only certain unperturbed states have 
strong admixtures into the given eigenstate. 

Such eigenstate structure is quite different from the 
case of conventional BRM. For a better understanding 
of these spiked fluctuations we study the local spectral 
density pw introduced by Wigner H and analyzed in 
BRM with linearly growing diagonal corresponding to 
conservative systems : 



Pw(E — E n ) 



A 



\i> x (n)\ 2 5(E - E X ) 



(2) 



The function pw characterizes the average probability 
P{\ijj\{n)\ 2 ) — pw{E — E n ) of eigenfunction ip x { n ) on 
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site n with energy E n = H n ^ L , where A is the eigenvalue 
index and n marks the original basis. Our numerical 
investigations in a wide range of parameters (20 < b < 
2000, 201 < N < 4001 and 1.5 < W b < 40) both in 
localized (l sb <C N) and delocalized (l sb ^> N) regimes 
show (see Fig. 2) that pw is well described by the well- 
known Breit-Wigner distribution pw = Pbw'- 

Pbw {e - E n ) = M{E _ + r2/4) ; r = ^ (3) 

where T is the distribution width. This distribution re- 
mains valid in localized and delocalized regimes under 
the condition that T is much less than the energy width 
SE f» 1 of H n . n i at Wb = 0. Usually the distribution 
Pbw appears in such physical systems as nuclei and com- 
plex atoms Q where due to energy conservation the di- 
agonal term rj n grows linearly with n that corresponds 
to a finite level density pe- In this case the width is 
r = 2-kpe < H^ n , > [||]. In our case for W b > 1 all 
eigenenergies are homogeneously distributed in the finite 
interval [— Wb, W>] and for full matrices with b = N/2 
we can use the above expression with pe — N/2W b 
which gives T in (3). For b <C N according to H one 
should replace pe by the density of directly coupled 
states p c = b/Wb that leads to the same expression for T. 
The theoretical formula for T, independent on b and N, 
is in a good agreement with our numerical data (Fig. 2). 
The independence of T on b and N makes our case quite 
different from the case of full matrix (1) studied before 
in§. 

For Wb 3> 1 the width of the Breit-Wigner peak is 
small and therefore according to (3) the probability on 
nearby levels is a strongly fluctuating spiked function. 
This spike structure of eigenfunctions can be character- 
ized by the IPR £a = (J2 n IV'a("-)| 4 ) _1 which counts the 
number of spikes independently on the distance between 
them. In the case of full matrices with b — N/2 the num- 
ber of spikes can be estimated as the number of states 
in the interval T that gives the average value of IPR 
£ =< 1/Ca > _1 ~ Pe^ ~ N/(2W b 2 ). The same estimate 
can be also used in the delocalized regime l sb 3> N with 
b <C N. Of course this estimate is valid only when the 
number of states in the width T is much larger than one, 
that implies £ > 1 or W b < V~N. 

The numerical results for the dependence of IPR on 
Wb in the delocalized regime are presented in Fig. 3. They 
demonstrate that for sufficiently large full matrices (N — 
4001) this dependence approaches to the above estimate. 
However the convergence is rather slow so that for smaller 
N values one has approximately £ ~ N/W b a where the 
exponent a slowly changes with N. For example a ~ 1.7 
for N = 2001. We attribute this very slow approach to 
the asymptotic value of a — 2 to the quite restricted 
range of W b variation. Indeed on one side the width of 
the Breit-Wigner peak should not exceed the width of the 
energy band for W b = that gives W b ^> 1. On the other 



side one should have W b <C Vn. Another restriction 
appears for band matrices with b < N/2 namely l sb ^> N. 
The data for this case (Fig. 3, full squares) show that 
for not very large W b the IPR is close to the regime of 
full matrices while for large W b one enters the localized 
regime l sb <C N which should be studied separately. 

It is interesting to note that in the delocalized regime 
even for W b 1 many levels are coupled by interaction 

if PB r«iv/(2Wb 2 )»i. 

Therefore, one would expect that for W b < W b cr »s 
(N/2) 1 ^ 2 the level spacing statistics P(s) will be the same 
as in the case of Gaussian orthogonal ensemble (GOE) 
ifjof . These expectations are not so evident since the 
spiked structure of eigenfunctions apparently should lead 
to a decrease of overlapping matrix elements between 
eigenfunctions. However, our numerical results for ma- 
trices with N < 8000 show that P(s) remains close to 
GOE for 1 < Wb < W b cr . They are also in agreement 
with the numerical results |ll[] for full matrices of smaller 
sizes showing that the transition border in W b between 
Poisson and GOE statistics scales as TV 1 / 2 . The question 
about other statistical properties of levels in the regime 
1 < Wb < W b cr remains open. 

For the localized regime in the above estimate of £ one 
should replace N by l sb since only levels in the interval of 
one localization length can contribute to the IPR. This 
gives the expression 

£ « l sb /2W b p - 2 « 6 2 /4Wf; (3 = 4 (4) 

which is valid for £ > 1 (1 < W b < Vb). The last 
condition together with l sb <C N gives strong restrictions 
for the numerical simulations (1 <C W b <C TV 1 / 4 ). 

Our results for this localized case are presented in 
Fig. 4. The data can be empirically fitted by £ ~ b 2 /W b 
with (3 ss 3 which differs from the theoretical value (3 = 4. 
We attribute this difference to the fact that we are not 
far enough in the asymptotic regime of large b and W b . 
Indeed, for W b > b 1 / 2 one enters in the perturbative 
regime and the deviations from a power law becomes ev- 
ident. We also checked that the probability distribution 
P(\if}\(n)\ 4 ) is proportional to p 2 BW that gives additional 
grounds for the theoretical power (3 = 4. However, the 
simulations with large enough values of parameters b, W b 
requires too large matrix sizes being beyond our numer- 
ical abilities. The numerically found value f3 > 2 implies 
that the number of peaks is smaller than the localization 
length l sb w 6 2 /(2W 2 ) which determines the asymptotic 
exponential decay of the eigenstates. It would be desir- 
able to have a more rigorous theoretical derivation of the 
IPR dependence on parameters in the localized regime. 

The above results show that the SBRM (1) has many 
features similar with the photonic localization in a molec- 
ular quasicontinuum Jl2] | as it was remarked in ||. Ac- 
cording to this analogy, the number of levels in one- 
photon transition (size) is of the order b and the den- 
sity of coupled states is 6/2Wb. However, in the photonic 
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model the levels are ordered in energy in a growing way 
that leads to a chain of equidistant Breit-Wigner peaks in 
an eigenstate @. For the SBRM (1) all levels are mixed 
in the energetic interval and the Breit-Wigner peak is 
hidden. 

Let us now discuss the consequences of the result (4) 
for the TIP model. According to the relation between 
the parameters of SBRM and TIP given above we obtain 
from (4) the expression for the IPR £ c in the TIP model: 

£ c ~ (U/V)% 2 > 1 (5) 

This result can be also derived directly from the den- 
sity of states inside the localization length interval l c 
(pE ~ hh/V) and the transition rate T c ~ U 2 /(Vli) ob- 
tained in ||. Indeed, the number of levels in the Breit- 
Wigner peak is T c pE ~ £ c that gives (5). This result 
shows that the number of noninteracting eigenstates £ c 
contributing in the eigenfuction is quite large for U ~ V 
while for (U/V) 2 li <C 1 this number is order of 1. How- 
ever, the value of £ c at U ~ V is much less than the 
number of unperturbed states AiV contributing to the 
TIP eigenfunction in the unperturbed lattice basis. This 
number determines the IPR ^ max « AN ss Z e Zi ~ if. 
The difference between £ c and £ max shows that the non- 
interacting eigenbasis represents the real eigenfunctions 
in a much better way. It also stresses the fact that the 
IPR value is not basis invariant. 

From the difference between £ c and ( mm it is possi- 
ble to conclude that the coherent propagation of TIP 
goes by rare jumps of size 1% between the states with 
approximately constant sum of noninteracting energies 
E s = e n + e n r. The distribution over E s should have 
the Breit-Wigner form with the width r c . The length 
of propagation by such jumps is l c ~ l\ 2 ^s> l\. Due to 
this hidden Breit-Wigner distribution the IPR £ c in the 
basis of noninteracting eigenstates is proportional to l\ 2 
instead of "naive" ^i 3 . For the case of TIP with M trans- 
verse channels one should replace l\ in (5) by Mix with 
li oc M being one-particle localization length. 

If one-particle motion is ergodic in a <i-dimensional sys- 
tem of size L < l\ then its eigenfunction contains about 
N\ ~ L d components. The matrix element of inter- 
action is then U s ~ U/Ni 3 ^ 2 ||, the density of cou- 
pled states p c ~ Ni 2 /V and the Breit-Wigner width 
r c ~ U 2 p c ~ U 2 /N{V for U < V is less than one- 
particle level spacing Ai ss V/N\. Thefore, it is possible 
that a concept of pairs formed by TIP can be useful even 
in the ergodic samples with L < 1%. In some sense, for 
A 2 < T c < Ai, where A 2 « V/N x 2 is two-particle level 
spacing, one can at first average over fast one-particle mo- 
tion and after that analyse the slow pair dynamics with 
typical time scale l/r c . In the ergodic regime L < li the 
IPR is £ c ~ T c p c ~ Ni{U/V) 2 < N 2 and according to 
the discussed above properties of P(s) in SBRM and the 
result the GOE statistics for TIP should be observed 



for £ c > 1. For L 3> h the strong enhancement of inter- 
action (p c T c ~ li d (U/V) 2 ^> 1) leads to derealization 
of the TIP pairs in d > 3 below one-particle Anderson 
transition when noninteracting particles are well local- 
ized || . 
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FIG. 1. Localized eigenfunction of a SBRM with Wb = 7, 
b = 100 and N — 2001. The solid line indicates the expo- 
nential localization with i s (, ~ 171 in agreement with results 
obtained in [3], eq.(3). 
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FIG. 2. Local spectral density which determines the aver- 
age probability on a given site P(\ip\(n)\ 2 ) — pw(E — E„) for 
b = 100, W b = 5, N = 201 (triangles, 20 realisations of disor- 
der) and N = 2001 (squares, 2 realisations of disorder). The 
solid line gives Breit-Wigner distribution (3) with V = 0.21. 
The inset shows the dependence of T on W b _1 : points are 
numerical data (N = 1001,6 = 100), straight line is theory 
(3). 



FIG. 3. IPR £ normalized with its limit value for the GOE 
case N/3 vs. Wb in the delocalized regime for N = 2001, 
b = 300 (full squares), and full matrices with N = 251 (open 
squares), N = 501 (open triangles), N = 1001 (x), N = 2001 
(full triangles) and N = 4001 (full circles). Dashed line shows 
the fit for full circles with a = 1.75 ± 0.03; solid lines shows 
theoretical slope a = 2. 



FIG. 4. Dependence of £/b 2 on Wb in localized regime : 
N — 2001, b — 50 (full squares) and b = 80 (open squares); 
N = 4001, b = 50 (full circles) and b = 100 (open cir- 
cles). Dashed line shows the slope from fit for open circles 
(/3 = 3.0±0.1) and solid line indicates theoretical slope f3 = 4. 
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